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We report an investigation of the structure of ^^C nucleus employing a newly developed 
configuration-mixing method. In the three-dimensional coordinate-space representation, we gen- 
erate a number of Slater determinants with various correlated structures using the imaginary-time 
algorithm. We then diagonalize a many-body Hamiltonian with the Skyrme interaction in the space 
spanned by the Slater determinants with parity and angular momentum projections. Our calcu- 
lation reasonably describes the ground and excited states of ^^C nucleus, both for shell-model-like 
and cluster-like states. The excitation energies and transition strengths of the ground-state rota- 
tional band are well reproduced. Negative parity excited states, IJ", 2^, and 3]~, are also reasonably 
described. The second and third 0"*" states, 0^ and Oj, appear at around 8.8 MeV and 15 MeV, 
respectively. The Oj state shows a structure consistent with former results of the a-cluster mod- 
els, however, the calculated radius of the Oj state is smaller than those calculations. The three-a 
linear-chain configuration dominates in the O3" state. 

PACS numbers: 21.10.Dr,27.20.-|-n 
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I. INTRODUCTION 

Light nuclei show a variety of structures in the ground 
and excited states characterized by different correlations. 
The nuclear shell model usually provides a reasonable de- 
scription for the ground and low-lying states. However, it 
is difficult for the shell model to describe well-developed 
cluster states which appear in excited states of light nu- 
clei. The appearance of the cluster structures is inti- 
mately related to the small separation energy of clusters. 
This is called the threshold rule and has been described 
by the so-called Ikeda diagram [1]. 

For theoretical descriptions of cluster states, micro- 
scopic and semi-microscopic cluster models have been 
extensively developed in the past [2-5]. The resonating 
group method (RGM) [6, 7] assumes a product form for 
the wave function composed of the internal wave func- 
tions of clusters and the inter-clustcr wave function, tak- 
ing fully account of the antisymmetrization. The gen- 
erator coordinate method (GCM) was also successfully 
applied to various cluster motions assuming harmonic os- 
cillator shell-model wave functions for clusters, as well 
[8]. In these microscopic models, the existence of clusters 
is assumed from the beginning. To understand mecha- 
nisms of emergence and disappearance of cluster struc- 
tures, one should start with a model which does not as- 
sume any existence of clusters. Studies with the antisym- 
metrized molecular dynamics (AMD) method [9-11] have 
contributed to substantial advances in this direction. In 
most calculations with microscopic cluster models and 
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AMD, effective nucleon-nucleon forces are used. In the 
Fermionic Molecular Dynamics (FMD) method which is 
closely related to the AMD, a more realistic force pro- 
duced by the unitary correlation operator method has 
been employed [12]. 

In last two decades, there have been significant ad- 
vances in theoretical descriptions of light nuclei start- 
ing with realistic nucleon-nucleon force, so-called the 
ab-initio approach. The Green Function Monte Carlo 
(GFMC) approach has been successful to describe the 
ground and the low-lying excited states of light nuclei 
[13]. The GFMC calculation describes the two a cluster 
structure of ^Be in the ground state [14]. The no-core 
shell model (NCSM) has also been successful for the de- 
scription of ground and low-lying excited states [15, 16]. 
However, descriptions of cluster states in the NCSM have 
not yet been satisfactory. Recently there are a number 
of attempts for the ab-initio description of cluster struc- 
tures in excited states. For example, a lattice calculation 
for the Hoyle state has been reported [17]. The no-core 
shell model combined with Monte-Carlo basis generation 
method has also been applied [18]. 

There are two important ingredients in the ab-initio 
descriptions of nuclear structures. One is to start with a 
Hamiltonian with a realistic nucleon-nucleon force that 
has a short-range repulsive core. The other is to obtain 
fully convergent solutions for the many-body Hamilto- 
nian. Since cluster structures are characterized by long- 
range spatial correlations, simultaneous descriptions of 
both long- and short-range correlations are required in 
the ab-initio calculations of cluster states. This makes 
the problem computationally very challenging. 

In this paper, we focus on the latter aspect of the 
above-mentioned problem, namely, on obtaining fully 



convergent solutions for a given many-body Hamiltonian, 
taking into account a variety of long-range correlations. 
We start not with a realistic nucleon-nucleon force but 
with an empirical effective interaction, the Skyrme force. 
We use a newly developed method which was reported 
previously [19]. We apply the method to ^^C nucleus, 
then, examine descriptions of cluster structures. Since 
the Skyrme interaction is determined so as to reproduce 
nuclear properties of wide mass region, our calculation 
contains no empirical parameters specific to ^^C. 

Among light nuclei, the ^^C nucleus is one of the most 
interesting systems for the reasons described below. In 
the jj coupling shell-model picture, the ground state 
wave function should be dominated by the P3/2 closed 
shell configuration. Indeed, the self-consistent Hartree- 
Fock solutions with most Skyrme interactions show a 
spherical shape with the P3/2 closed shell configuration. 
However, the ^^C nucleus is known to show a rotational 
band structure built on the ground state, indicating a 
deformed intrinsic shape in the ground state. In excited 
states, a variety of cluster structures are known to ap- 
pear. The Oj state just above the three-a decay thresh- 
old is an important resonant state for the triple-a fu- 
sion reaction, and is known as the Hoyle state [20]. It 
has been found recently that this state is well described 
by a Bose condensed wave function of three-a particles 
which is called the THSR wave function [21, 22]. The 
appearance of three-a linear-chain structure in excited 
states was suggested by Morinaga in 1966 [23]. Recent 
microscopic cluster models predict that the 0^ state is a 
candidate for the linear-chain like structure [11, 12]. 

This paper is organized as follows. In section II, we 
explain our method. In section III, we first show that 
we obtain convergent results for some low-lying states 
for the Skyrme Hamiltonian. We then describe in detail 
the calculated results for ^•^C. In section IV, we compare 
our results with conventional microscopic cluster model 
calculations by introducing cluster wave functions in our 
calculation. In section V, a summary is presented. 



II. FORMULATION 

In this section, we present our formalism which was 
proposed in Ref. [19]. It is composed of three steps: We 
first generate a number of Slater determinants (SDs) in 
a stochastic way. These SDs are expected to span a suffi- 
ciently large model space to describe excited states with 
various cluster structures as well as low-lying states with 
shell-model-like structures. We then perform parity and 
angular momentum projections for the SDs. Finally we 
superpose them to diagonalize the Skyrme Hamiltonian. 
Below we describe these three steps of our formalism in 
order. 

Before presenting our formalism, we briefiy describe 
numerical aspects in the present method. For the en- 
ergy functional, we employ the SLy4 parameter set of 
the Skyrme interaction unless otherwise specified. To 



describe single- particle orbitals of the SDs, we employ a 
representation of the three-dimensional (3D) Cartesian 
grid. This representation allows us a flexible descrip- 
tion of single-particle orbitals in arbitrary nuclear shapes. 
The grid spacing is taken to be Ax = Ay = Az = 0.8 
fm, and all the grid points inside a sphere of radius, 
^max ~ 8.0 fm, are adopted. 



A. Generation of Slater determinants 

The first step is to generate and select a sufficient num- 
ber of SDs, {$i} (i = 1,---,M), which span a large 
model space to describe various kinds of long-range corre- 
lations. To this end, we make use of the imaginary-time 
method with the Skyrme interaction. The imaginary- 
time method is usually used to obtain a ground-state 
solution in the mean-field calculation. Here, we use it for 
generation of many kinds of collective surfaces. 

We start with initial SDs whose single-particle orbitals 
are described by the Gaussian wave packets. 
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The centers of the Gaussian wave packets, Ri{a), are set 
by random numbers generated under the condition. 



\R^{(T)\ <i?„ 



1 fm. 



(2) 



The width parameter a is taken to be 2 fm. 

Then, we start the imaginary-time iterations with 
these initial SDs. During the imaginary-time iterations, 
we set constraints to place the center of mass at the ori- 
gin and to make the principal axes of nuclei parallel to 
the Cartesian axes. After the sufficient number of itera- 
tions, it reaches the self-consistent ground-state solution, 
namely Hartree-Fock (HF) state. However, before reach- 
ing the HF solution, there appear a number of configu- 
rations which show various cluster structures and other 
configurations important for low-energy nuclear dynam- 
ics. Thus, we pick up and store these SDs which will 
be used for the configuration mixing calculation. We re- 
peatedly perform the imaginary-time iterations starting 
with different initial SDs to obtain a sufficient number of 
SDs (typically 50) of many kinds of correlations. 

In order to span a wide model space by the limited 
number of SDs, the actual selection of the SDs is achieved 
as follows: The first SD adopted in the basis set {$i} 
is the HF state, $1 = 'I'hf- The second and following 
SDs are generated as follows: During the imaginary-time 
iteration, the energy expectation value decreases mono- 
tonically. We examine the energy expectation value in 
a regular interval, typically every 100 iterations. We do 
not adopt the SD until the energy expectation value is 
less than 30 MeV above the HF ground-state energy. Af- 
ter the energy expectation value becomes less than this 
threshold, we examine the overlap between a current SD 
and previously selected ones. We denote the SDs which 




23th 



34th 



39th 



50th 




FIG. 1. Contour plots of nuclear densities of the stored SDs for 
the panel. Units of vertical and horizontal axes are fm. 



C. A sequential number of the SD is indicated in the top of 



have been already selected as $i(i = 1, • • • ,N), and de- 
note the current SD as $c- We calculate the quantities 



|($,|F"i?„|$c)| 



($,|P-|$,)|J|($,|P-|$,)| 



ii = l,---,N), (3) 



shows a ^Be+a like structure. We thus observe that the 
present procedure efficiently produces SDs with various 
cluster structures in an automatic manner. 



B. Projections of parity and angular momentum 



where P^ is the parity projection operator and Rn(n ~ 
1 , ■ • • , 24) are operators of rotations and inversions which 
may be easily achieved by the changes of the coordinate 
axes. If the maximum value of the above overlap values 
between $c and $i(i = 1,---,A^) is less than 0.7, we 
adopt the SD $c as a new member, $7v+i = $c. When 
the imaginary-time iteration reaches the self-consistent 
solution, we generate a new initial SD of the form of 
Eq. (1) and again start the imaginary-time iterations. 
We then repeat the procedure to select a new SD. We 
repeat the procedure until a sufficient number of SDs is 
obtained. Typically a few SDs are adopted during one 
imaginary-time iterations. As the number of stored SDs 
increases, it is more and more difficult to find the new 
SD which satisfies the overlap condition. In this way, we 
store SDs which include important correlation effects and 
which are linearly independent to each other as much as 
possible. 

We show density distributions of several SDs gener- 
ated in this procedure in Fig. 1. The numbers assigned 
to the figures, 1, 3, 14, 21, 23, 34, 39, and 50, simply in- 
dicates the adopted ordering. The first one which shows 
a spherical shape is the HF solution for the ground state. 
Other SDs in Fig. 1 show a variety of cluster structures. 
For example, $3 shows an equilateral triangular three- 
a structure, $14 shows a thrce-a linear-chain, and $23 



The SDs prepared by the method in Sec. II A are, in 
general, not eigenstates of parity and angular momen- 
tum. To calculate matrix elements between eigenstates 
of parity and angular momentum, wc apply the projec- 
tion method. The projection operators are given as usual 
by 



p- = l(i + (-irA), 



(4) 
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where Pr is the space inversion operator, i?(a, /3, 7) is the 
rotation operator for the Euler angles, a, /3, and 7, and 
^iiK i^ the Wigner's Z3-function defined by 



i?(a,/3,7) 



^-iaJ^^-iPJy^-i-fJ^ 



^-iaM jj 



'■MK 



(fi)e 



-i^K 



(6) 

(7) 



where J, M, and K are the total angular momentum, its 
projection onto the laboratory z-axis, and its projection 
onto the body-fixed z-axis, respectively. 

Wc define the norm and Hamiltonian matrix elements 
between the projected SDs |<f>i) and 1$^) as 
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Here, we use the formula, 
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where e;^'^ and v^'^^ are eigenvalues and eigenvectors of 
the norm matrix, n'l^^j^,. We then construct a space 
spanned by the eigenvectors with the eigenvalue efj > 
10~2, and define the normalized basis functions 
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In Eqs. (8) and (9), we need the rotation of wave func- 
tions. It is achieved by successive operations of small- 
angle rotations. For example, the rotation of a wave 
function (p over an angle 7 around the z-axis is achieved 
by successive rotations of a small angle, A7 = 7/^7, ^7 
times. 



e-'^^'\(l,) = (e- 
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To achieve the small angle rotation, we employ the Taylor 
expansion method. 
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Typically, we take -/V,„ax — 4 and A7 = ^. 

The integrals over Euler angles in Eqs. (8) and (9) are 
evaluated as follows: Those over a and 7 are achieved 
by the trapezoidal rule, taking 18 uniform grid points for 
[0, 27r]. The integral over /3 is achieved with the Gauss- 
Legcndrc quadrature, taking 30 grid points for [0,7r]. 



C. Configuration mixing 

The final procedure is diagonalization of the many- 
body Hamiltonian in the space spanned by the selected 
SDs. In Sec. II A, the SDs have been screened by their 
linear independence. However, calculating eigenvalues of 
the norm matrix for the SDs after the parity and angu- 
lar momentum projections, we find a number of eigen- 
values very close to zero or even slightly negative. The 
norm matrix is positive definite by definition. However, 
since we make numerical approximations in evaluating 
the norm matrix, it could contain negative eigenvalues. 
The approximations include use of the formula for the 
product of the projection operators, Eq. (10), which is 
no longer exact if the integral over Euler angles is eval- 
uated by the numerical quadrature. We also employ the 
3D Cartesian grid representation for the orbitals in which 
the rotational symmetry holds only approximately. 

Inclusion of those configurations of very small norm 
eigenvalues would lead to numerical difficulties. There- 
fore, we reduce the number of configurations according to 
the following procedures. First, we perform diagonaliza- 
tion in the (2 J + l)-multiplet with different K quantum 
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After achieving the above procedure for all SDs, we de- 
fine the following norm matrix between basis functions 
belonging to the different SDs, 
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Using this matrix, we examine the linear independence 
of the basis functions and reduce the number of basis as 
follows. 

1. Calculate eigenvalues of 2 x 2 matrices composed 
of every possible pair of basis functions, (i\V\) and 
(i2V2)- If the smaller eigenvalue is less than 10~^, 
we remove the basis function with a smaller e/j^. 

2. Calculate the eigenvalues of 3 x 3 matrices com- 
posed of three basis functions (iii^i), [12^2) and 
{i-iV-^. If the smallest eigenvalue is less than 10~^, 
we remove one of the three basis functions in the 
following procedure. We calculate eigenvalues of 
three 2x2 submatrices composed of all possible 
pairs of these three states, to find the pair whose 
smaller eigenvalue is the largest among the three. 
Then we remove one of the basis functions of (iiJ^i), 
(12^2) and {i-iV-i) which does not belong to that 
pair. We repeat the procedure for all possible com- 
binations of three basis functions. 

3. Finally we calculate eigenvalues of the norm ma- 
trix fi'l'^-j with basis functions which survived in 
the previous two screening steps. If we find the 
eigenvalue smaller than 10~'^, we remove one ba- 
sis function in the following way. Denoting the 
number of basis functions as N , we construct the 
{N — 1) X {N — 1) submatrices removing one ba- 
sis function $/j^ from the N basis. Apparently, N 
different choices of {iv) are possible. We then calcu- 
late the smallest eigenvalue of the {N — 1) x {N — 1) 

submatrix, A^^^j^. Among X^^^l with different (w), 

we find the largest one, A^!,;|^ , and remove the ba- 
sis function ^fi^^i ■ In this way, the number of basis 
{$fj^ } is reduce by one, from iV to iV — 1. We 
repeat this process until the smallest eigenvalue of 
the norm matrix becomes larger than 10~^. 

After removing the overcomplete basis functions in this 
procedure, we achieve the configuration mixing calcula- 



tion. Denoting the nth energy eigenstate as 

l^;f) = E/-""l*-), (16) 

the generahzed eigenvalue equation for the energy eigen- 
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III. RESULTS 
A. Convergence of results: Statistical treatment 

In principle, the configuration mixing calculation with 
a sufficient number of SDs should provide unique and 
convergent energy levels. However, as we described in 
Sec. II, superposing a large number of non-orthogonal 
SDs causes numerical difficulties. In the present calcu- 
lations, we adopt 50 SDs for the configuration mixing 
calculation. It is difficult to increase this number. Fur- 
ther increase of the number of SDs may produce unphys- 
ical solutions whose energies are a few tens of MeV lower 
than the ground state energy of the HF solution. This is 
possibly due to accumulations of numerical errors by the 
violation of rotational symmetries in the 3D grid repre- 
sentation, insufficient accuracy in numerical quadratures, 
and so on. 

Because of the difficulty, we will not attempt to exam- 
ine the convergence of the energy levels by increasing the 
number of SDs. Instead, we prepare several sets of the 
SDs and calculate energy levels for each set. If the cal- 
culated energy levels are close to each other among the 
different sets of SDs, one may conclude that the calcu- 
lated energy levels are reliable. In practice, we prepare 
ten sets, each of which is composed of 50 SDs. The ten 
sets of SDs are prepared in the procedure explained in 
Sec. II A. Different seeds for the random numbers, which 
are used to prepare initial states in Eq. (1), are employed 
to generate the different sets. 

In Fig. 2, we show the energy levels of ^^C nu- 
cleus for the ten sets calculated in the procedure ex- 
plained in Sec. II. The energy levels are shown for J^ = 
0+,l+,2+,3+,4+,l-,2- and 3". 

Let us first examine calculated energy levels of J'^ = 
0+. The lowest level is located around -95 McV. The 
difference of energies among the ten sets is smaller than 1 
MeV. The second excited state appears around -86 MeV. 
The difference among the ten sets is again around 1 MeV. 
The third excited state appears around -81 MeV. We 
may state that the energies of these three lowest states 
arc calculated reliably, since the variation is rather small. 
However, energies of fourth excited state do not show a 



good convergence. For example, the energy levels of 2nd 
set give the energy at around -79 MeV, close to the 3rd 
state. However, the energy in the 7th set is substantially 
high, approximately -76 MeV. We thus conclude that we 
can obtain reliable excitation energies and wave functions 
for the lowest three levels for J'" = 0+. 

The energy levels of J^ = 2+ in Fig. 2 indicate that the 
energies of the lowest four states are reliable with a small 
variation. For J'^ ~ 3+ and 4+ states, the lowest two 
states may be reliable. However, the calculated energies 
of J'^ = 1+ levels show strong variation among the ten 
sets even for the lowest level. This may be due to the 
fact that the J'^ = l"*" components of the wave function 
disappear in early stages of the imaginary-time iterations, 
since components of high-lying levels quickly decay by 
the imaginary-time propagation. For the negative-parity 
levels, only the lowest level for each J'^ may be reliable. 
The energies of second lowest levels show a large variation 
among the ten sets for J'^ = 1~, 2^ and 3^. 

For physical quantities such as energies, transition 
strengths, and radii, we calculate statistical averages and 
standard deviations among the ten sets. The average en- 
ergy for the n-th level of J'^ state is defined by 
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where i specifics a set among the ten sets, Ng = 10. The 

average excitation energies are calculated as E^ — Eq , 
which will be shown in Figs. 3 and 4. We also calculate 
standard deviation of the energies which will be shown 
by the error bar in the figure. The standard deviation is 
defined by 
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The average values and the standard deviations for the 
transition strengths and radii are evaluated in the same 

way. 



B. Energy levels 

We show calculated excitation spectra of even and odd 
parities in Figs. 3 and 4, respectively. In the figures, 
energies averaged over ten sets are shown with error bars 
as the standard deviation. Our calculated results are 
compared with measurements [27-31] and other theories, 
AMD [11], GCM [24], RGM [25], and NCSM [26]. 

In the Skyrme-HF calculation, the binding energy of 
^^C is 90.6 MeV, in reasonable agreement with the mea- 
sured value, 92.2 MeV. In our configuration mixing cal- 
culation, the correlation energy is 4.7 ± 0.2 MeV. The 
ground-state energy including the correlation is 95.3 ± 
0.2MeV, slightly lower than the measured value. 

Calculated excitation energies of the ground-state ro- 
tational band arc in good agreement with measurements. 
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FIG. 2. Energy levels of ^^C nucleus for J'^ = 0^, 1^, 2"*", 3^, 4"'", 1 ,2 and 3 . Calculations employing ten different sets of 
SDs are shown. See the text for details. 
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FIG. 3. Excitation energies of positive parity are shown. The 
energies are obtained by averaging over ten configurations. 
The standard deviations of the energies are also shown by 
error bars. We also show the results of AMD [11], GCM [24], 
RGM [251 and NCSM [261. See the text for details. 



In the nice reproduction of the rotational energy lev- 
els, the configuration mixing is essential since the ground 
state solution in the HF calculation is spherical for ^^C 
with the SLy4 interaction. The excitation energy of 
2\ state is well reproduced by the present calculation, 
AMD, and NCSM. However, microscopic a cluster mod- 
els (GCM and RGM) provide too low excitation energies. 
The former models (present, AMD, and NCSM) take into 



FIG. 4. Excitation energies of negative parity states. The 
same explanation as that in Fig. 3 applies to others. 



account the spin-orbit interaction, while it is not included 
in the latter (GCM and RGM) in which existence of the 
three a clusters are assumed. This suggests that a proper 
inclusion of the spin-orbit interaction is important for the 
good description of the ground rotational band. 

The Oj state, known as the Hoyle state, has been at- 
tracting much attention recently since it has been shown 
that this state may be understood as the Bose condensed 
state of three a particles [22]. Our calculation gives a 
reasonable description for this state, although the exci- 
tation energy is slightly overestimated by about 1 MeV. 



Properties of this state will be discussed in the following 
subsections. Although recent ab-initio approaches have 
been successful for the ground-state rotational band, a 
satisfactory description for the 0^ state has not yet been 
made. For example, the NCSM cannot describe this 
state adequately [26] . Recently, attempts of ab-initio de- 
scription for this state have been undertaken by several 
groups. A nuclear lattice calculation for this state has 
been reported in Rcf. [17]. 

Recently, a new 1^ state has been found experimen- 
tally at 9.6 ± 0.1 McV excitation energy with a width 
of 0.6 ±0.1 MeV [28, 32] This state was interpreted as 
the excited state built on the Oj state. In our calcula- 
tion, two 2+ states, 2 J and 2^, appear just above the 
Oj state. However, as we discuss in Sec. IIIC, these two 
states, 2 J and 2j^, in the present calculation, seem not 
to correspond to rotationally excited state of the Hoyle 
state. 

In Fig. 3, three states, 0;^, 2^, and 4^, follow a rota- 
tional energy sequence. Small standard deviations of the 
energies of these states indicate the reliability of the cal- 
culation. As will be discussed in Sec. IIIC, these states 
are connected by strong B{E2) transitions. In Sec. IIIF, 
we will show that this band corresponds to a three-a 
linear-chain state. 

For the negative parity states, we have obtained solid 
results only for the lowest energy state for each J^ sector 
(Sec. Ill A). Our calculation reproduces the measured 
order of the three states, 3j~, Ij^ and 2^ . However, the 
excitation energies are higher than the measurements by 
2-3 MeV. 

We would like to stress that our calculation includes no 
empirical parameter specific to the system, ^^C nucleus. 
We employ the SLy4 parameter set which is determined 
to reproduce nuclear properties of whole mass region. 
This is in contrast to cluster model calculations where 
nuclear force parameters are often adjusted for respective 
systems. We also do not employ any effective charges in 
evaluating the transition matrix elements shown below. 

Finally we mention how the calculated energy levels 
depend on the choice of the Skyrme interaction. In Fig. 5, 
we show the excitation energies of positive parity states 
with different parameter sets of the Skyrme interaction, 
SLy4, SkM*, and SHI. The same set of SDs (No. 1 in 
Fig. 2) is employed in all calculations. The correlation 
energies in the ground state are shown as well inside the 
parentheses. The comparison shows that basic features 
of the spectra do not depend much on the choice of the 
Skyrme parameters. For example, the ground rotational 
band is described reasonably by all three parameter sets. 
The position of Oj state does not change much. There 
appear rotational band in three calculations starting with 
0;^ state. We thus conclude that the excitation energies 
are not sensitive to the choice of the Skyrme parameters 
for almost all the states below 15 MeV. 
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FIG. 5. Energy levels of positive parity employing different 
parameter sets of Skyrme force, SLy4, SkM* and SIII. The 
number in parentheses is the correlation energy in the ground 
state, Ehf — Egs, in unit of MeV. In the calculation, the same 
set of SDs is employed. 



C. Transition strength 

Calculated B{E2), B{E3), and M{EO) values, the av- 
erage values and the standard deviations, are shown in 
Table I. In our calculated values, we do not employ any 
effective charges. The B{E2) transition strength between 
2^ and 0^ states is well reproduced by our calculation. 
It is also consistent with results of other theories. The 
standard deviation is small, about 3%, indicating that 
our calculated value is well converged. 

The B{E2) transition between 0^ and 2^1 is calculated 
as 13.6 ± 1.2e^fm , which is in excellent agreement with 
the measured value, 13 ± 2e^fm . Other theories fail to 
reproduce the rate. In Ref. [11], it is argued that this 
transition strength is sensitive to the a-breaking effect. A 
good reproduction of this transition strength by our cal- 
culation indicates that our calculation reasonably takes 
account of the a-cluster components in the states. 

As mentioned in Sec. IIIB, there appear two 2+ states, 
2^ and 2|^, above the 0^ state in our calculation. These 
states might correspond to the 2"'" state at 9.6 MeV which 
was discovered recently [28, 32]. It was suggested to be 
a candidate of rotationally excited state built on the Oj 
state. In the present calculation, however, the B{E2) 
transition strength between 0^ and 2^ states is small as 
seen in Table I. The rate between 0^ and 2;^ states is 
also not very large. The B{E2) rate from Oj state is 
the largest for 2^ state which is regarded as rotationally 
excited state built on the 0;^ state, as will be mentioned 
below. These observations suggest that the 2^ and 2j^ 
states in the present calculation do not correspond to a 
rotationally excited state on the Oj state. 

As we discussed in Sec. IIIB, the states of 0;^, 2^, 
and 4^ follow the rotational energy sequence. The cal- 
culated transition strengths of B{E2; 2^ -^ 0'^) and 
B{E2;A2 — > 2^) arc very large. These results strongly 
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TABLE I. B{E2), B{E3) and M{EO) values of ^^C in units of e^fm*, e^fm'' and efm^ respectively. Experimental and calculated 
values are shown in the first and second column, respectively. For comparison, results of the AMD[11], GCM[24], RGM[25], 
NCSM[26], and THSR are shown. Values in THSR are calculated with the same model as Ref. [22]. The values indicated by * 



correspond to B{E2; 22 



05 



B{E2; 2t -> O;^) and B{E2;4+ -> 2^) in Ref. [11]. See text for details. 



support that these states indeed constitute a rotational 
band. In Sec. IIIF, we will see that this band is domi- 
nated by the threc-a hnear-chain structure. 

In the AMD calculation [11], intense B{E2) values are 
reported in the transitions among 0;^, 2^, and 4^ states. 
Since the states corresponding to 22 and 23 in our cal- 
culation seem not to be present in the AMD calcula- 
tion, we put these B{E2) values by AMD at the places 
of B{E2;2J -^ 0^) and B{E2;A+ -^ 2j) in Table I. 
In Ref. [11], these states are considered as the threc-a 
linear-chain states. The large B{E2) values are quali- 
tatively consistent with our result, though the absolute 
magnitudes of the transition strengths are much smaller 
in the present calculation. 

For negative-parity states, experimental data for 
B{E3;3^ -^ 0^) are available. The present calculation 
gives 77 ± 4e^fm , which is slightly smaller than the mea- 
sured value, 107 ± 14e^fm^. 

Finally, we discuss the M{EO) transition strength be- 
tween 0^ and O^j*" states. In the studies by cluster models, 
it has been argued that the magnitude of this transition 
strength reflects the spatial extension of the 0^ state [33] . 
Our calculated value, 4.5 ± 0.2efm , is slightly smaller 
than the measured value, 5.4 ± 0.2efm^. In contrast, 
microscopic cluster models and AMD have reported an 
opposite feature, slightly larger values, 6.6 — 6.7 e fm^, 
than measurement [11, 24, 25]. Experimentally measured 
value, 5.4 ± 0.2efm^ [27], is located between our result 
and those of the other calculations. 



D. Radii 

We next examine root-mean-square (rms) radii of the 
ground and excited states. Since our wave function does 
not allow an exact separation of the center-of-mass mo- 
tion from the internal one, we estimate an approximate 
correction for the radius due to the center-of-mass mo- 



tion, and subtract it from the calculated values. We as- 
sume a harmonic oscillator motion for the center-of-mass 
with the oscillator constant given by hoj = 41A~^^^ = 
17.9 MeV. The value of the correction in this model 
is estimated to be 0.07 fm in the harmonic oscillator 
shell model. The calculated radii after the correction 
are shown in Table II. 

Our calculated value in the ground state is 2.52 ± 0.01 
fm. This value is somewhat larger than the measured 
value, 2.31 ±0.02 fm. In the HF calculation, the radius is 
given by 2.44 fm. Our configuration mixing calculation, 
therefore, slightly increases the radius. Comparing with 
other theories, our value is larger than those of GCM and 
FMD, and is comparable to the value of AMD. 

For the 2^ state, our calculated radius is slightly larger 
than that of the ground state. Other theories report al- 
most the same or slightly larger radius for this state. 

For the Oj state, we find a significant difference be- 
tween the present calculation and the others. Our cal- 
culated radius is 2.73 ± 0.02 fm, which is larger than 
the radius in the ground state. However, this is much 
smaller than the other calculations which give more than 
3 fm [11, 22, 24, 25, 34]. In the recent AMD±GCM cal- 
culation [35], the radius of 2.9 fm was reported, similar 
to ours. It has been found that the radius of the 0^ 
state is quite sensitive to the spin-orbit interaction used 
in the AMD calculation [36] . The radius of the J state 
decreases as the strength of the spin-orbit interaction in- 
creases. This dependence is understood as follows [36]. 
If the strength of the spin-orbit interaction is weak, the 
ground state wave function contains substantial amount 
of the three-a-cluster component. Then, the 0^ wave 
function, which is dominated by dilute three-a compo- 
nents, spatially expands to ensure the orthogonalization 
to the ground state. As the spin-orbit interaction in- 
creases, the three-a component decreases in the ground 
state, which allows O2 wave function to include more 
compact thrce-a structure. This change results in de- 
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TABLE II. Mass rrns radii of the ground and excited states of 
^^C. The experimental data is taken from Ref. [38]. For com- 
parison, we show the resuhs of AMD [11], FMD [34], GCM 
[24] , RGM [25] and THSR. 



crease of the radius in Oj state. This mechanism may 
explain the discrepancy in the 0^ radius between our cal- 
culation and other theories. It should be noted again that 
our calculated value for the M{EQ) transition strength is 
smaller than the calculated values by other theories. Wc 
also note that an indirect measurement of radius for the 
0^ state using diffraction inelastic scattering [37] was re- 
ported, giving 2.89 ± 0.04 fm. 

For the Oj]" state, our calculated radius is 3.20 ± 0.05 
fm, which is much larger than the radii of 0^ and Oj 
states. This is again smaller than those by other models 
listed in Table. II, while it is similar to the value (3.26 fm) 
in Ref. [35]. 



E. Charge form factors 

A charge form factor from the initial state I J^, M^) to 
the final state | J/, M/) is defined as follows. 



\F.j,^jM^)?= 



1 1 

Z 2 J, + 1 



MiMj k 



(21) 



where Z is the proton number and q is the transferred 
momentum. Fp{q^) is a correction factor for the pro- 
ton size for which we employ Fp{q^) = exp(— 0^9^/6) 
with Op = 0.831 fm. To correct the center-of-mass mo- 
tion, we simply assume that the center-of-mass motion is 
separated and described by the harmonic oscillator wave 
function of the same oscillator constant, Tuo = 41A~^/'^ 
MeV, for both initial and final states. Thus, this leads 
to F^^{q'^) = exp(g262/2A) with 6 = 1.66 fm. 

In Fig. 6, we show charge form factors for the elastic 
(left) and inelastic 0^ -^ Oj (right) processes. Red solid 
curves show our results, blue dashed curves show the re- 
sults of AMD calculation [43] , and crosses with error bars 
show experimental results [39-42]. For the elastic form 
factor, we also show that of Skyrme HF solution in the 
ground state. 

In the small momentum transfer region q^ < 2 fm^^, 
the elastic form factor is well reproduced by the calcu- 
lation. For q^ > 2 fm~^, our calculation underestimates 



the form factor, though position of the dip at around 
3 fm~^ is reproduced well. The inelastic form factor 
for 0^ -^ Oj transition is underestimated for a whole 
momentum transfer region. The position of the dip at 
around 4 fm~^ is reproduced well. 

We show results by the AMD calculation. They are in 
better agreement with measured values, although the dip 
position in the elastic form factor is located at somewhat 
smaller q^ value. Microscopic cluster calculations also 
reproduce the form factors well [24, 25]. 

The underestimation of the elastic form factor at large 
q^ value indicates that the density in our calculation 
lacks high momentum component. Since the HF solu- 
tion gives a better description for the form factor at high 
momentum, the superposition of a number of Slater de- 
terminants turns out to increase the diffuseness in the 
nuclear surface, making the density distribution function 
p{r) smoother. As for the underestimation in the inelas- 
tic form factor of 0^ -^ J transition, a possible reason is 
the difference in the character of the wave functions be- 
tween the ground and Oj states. As we discussed in the 
radii, a rather small radius of O2 state in our calculation 
may indicate a small fraction of three-alpha component 
in the ground state. The inelastic form factor may be re- 
duced if the correlation structures are different between 
two states. It has been argued that the magnitude of this 
form factor at small q^ is quite sensitive to the radius of 
the Oj state [44]: the magnitude of the form factor at 
small q^ reduces as the radius of the 0^ state increases. 
Our result here is opposite, however. The radius of O2 
state in our calculation is smaller than those by cluster 
models, and the magnitude of the inelastic form factor is 
also small. 



F. Analysis for wave functions 

In order to clarify what kinds of correlations are in- 
cluded in the wave function after configuration mixing, 
^n^ , we calculate the overlap between the energy eigen- 
state and the projected single SD state, 
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(22) 



and find the SDs which have large overlap values with 
^'1^ ■ We show density distributions of the SDs to visu- 
alize the correlations included. 

In the following, we use the sequential number of the 
SDs which we assigned in Sec. II A, using the result of the 
first set of SDs in Fig. 2. We also show the K quantum 
number of the SD and the value of the overlap, P;/'^'*^, 
in Eq. (22). 
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FIG. 6. Squared elastic form factor for the ground state (left) and Of — >■ Oj inelastic form factor (right) are shown. Our 
calculated results are compared with the HF calculation with a single SD, measurements [39-42] and the AMD calculation [43]. 
Here, the result of the first set of SDs in Fig. 2 is used. 
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TABLE III. The sequential number, the K-vaiue, and the 
squared overlap value are shown for the SDs which dominate 
in the wave function of the ground rotational band. 
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FIG. 7. Contour plots of the density distributions of the 15th 
and the 4th SDs, which are the major components of the 
ground state rotational band. 



In Table III, we show the sequential number of the SDs 
which have large overlap values with the wave function 
of the ground rotational band, Of , 2'^, and 4^. The 
overlap values P^'^'*^ defined by Eq. (22) and K values 
are shown as well. Since the SDs are non-orthogonal, the 
sum of the overlap values is not equal to but much larger 
than unity. 

In the ground state Oj*", the 15th SD has the largest 
overlap, showing 90% for the overlap value. In 2f and 
4^ states, the 4th SD is the largest component and the 
15th SD is the second largest. To illustrate the nuclear 
shape of these two SDs, we show contour plots of density 
distributions of the SDs in the j/z, zx, and xy planes in 



Fig. 7. As seen from the figure, they both show oblate 
deformed shapes. 

The self-consistent HF solution is assigned to the first 
SD (number 1). We should note that it does not ap- 
pear in the top ten components of the ground state. Its 
overlap value with 0^ is about 70%. For ^^C, the self- 
consistent HF solution with the SLy4 interaction is spher- 
ical with P3/2 closed shell configuration. The spherical 
solution cannot describe the rotational band observed in 
the measurement. As shown in Fig. 3 and Table I, our 
calculation accurately reproduces the energy levels and 
the B{E2) transitions among the states of the ground 
rotational band. This good reproduction is achieved by 
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TABLE IV. The sequential number, the K-va.\ue, and the 
squared overlap value are shown for the SDs which dominate 
in the wave function of the negative parity states 1^, 2^, 3^ 
and 4j~. 



a superposition of SDs of deformed shapes. 




2. Negative-parity states 

In Table IV, we show sequential numbers of the SD 
which have large overlap values with the wave function 
of the negative parity states, 3j", Ij", and 2j". We find 
the 4th SD, which appears in the ground rotational band, 
also dominates in the negative parity states. Other SDs 
which dominate in the negative-parity states are 21, 29, 
and 47. 

We show density distributions of these three SDs in 
Fig. 8. All of these SDs have similar oblate shapes with 
three-a-likc structure. Close look at the densities reveals 
that the 4th and the 29th SDs have a compact configu- 
ration, while the 21tli and the 47th SDs show spatially 
more extended three-a configurations forming an obtuse- 
angled triangle. 



3. 0+, 2+ 

In Tabic V, we show sequential numbers of the SDs 
which have large overlap values with the wave function 
of the states, 0^, 2 J and 2^^. 

Wc first examine the Hoylc state, Oj. Compared with 
the cases of the ground rotational band and the negative- 
parity states, the maximum value of the overlap is rather 
small, less than 50%. This indicates that the superpo- 
sition of a large number of SDs is essential to describe 
the Oj state. This is consistent with the cluster-model 
calculations [24, 25] and the picture of the a condensed 
state for the Oj state [22]. 

We show in Fig. 9 the density distributions of the SDs 
which have the largest and next largest overlaps with 
the Oj state, namely, the 9th and the 28th SDs. These 



FIG. 8. Contour plots of the density distributions of the 21th, 
29th, and 47th SDs, which are major components of the neg- 
ative parity states, lj~, 2^, 3]~ and 4^. 
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TABLE V. The sequential number, the A'-value, and the 
squared overlap value are shown for the SDs which dominate 
in the wave function of the Ojm 2 J and 2^. 



SDs have a well developed cluster structures of three a- 
particles. 

Regarding the 2 J and the 2J^ states, wc find that a 
number of configurations contribute to these states, as in 
the case of 0^ state. The SDs in 2^ and 2;^ are more or 
less similar. However, they are very different from those 
in the Oj state. This is consistent with our observation 
that the B{E2) transition strengths between 0^ and 2^ 
states, and between Oj and 2^ states are rather small 



(see Sec. IIIC). 
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FIG. 9. Contour plots of the density distributions of the 
9tli and the 28th SDs, which are major components in the Oj 
state. 



FIG. 10. Contour plots of the density distributions of the 
30th and the 40th SDs, which are the major components of 



the O: 



and 4, states. 
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TABLE VI. The sequential number, the A'-value, and the 
squared overlap value are shown for the SDs which dominate 
in the wave function of the O^^m 2^ and 4 J.. 



4. Linear-chain states 

As seen in Sec. IIIC, 0;^, 2^^ and 4 J states are con- 
nected by the intense B{E2) values. In Table VI, the 
sequential numbers of the SD which have large overlap 
values with the wave function of these states are shown. 
The 30th, the 40th, and the 19th SDs are commonly in- 
cluded in the three states. We show in Fig. 10 the density 
distributions of the 30th and the 40th SDs. They clearly 
show a bended linear-chain structure of three a particles. 



IV. MIXING OF THREE-ALPHA 
CONFIGURATIONS 

Some of the present results in Sec. Ill are found to 
be qualitatively different from those of AMD and micro- 



scopic cluster models. For example, the radius of the 0^ 
state is much smaller in our calculation. The charge form 
factor at large momentum transfer is described much bet- 
ter by other theories than ours. These facts may indicate 
that the imaginary-time propagation may not sufficiently 
produce a certain class of a-clustcr wave functions. In or- 
der to check whether explicit inclusion of a-cluster con- 
figurations bring large changes in the current results, we 
perform configuration-mixing calculations including the 
wave functions similar to those employed in the micro- 
scopic cluster model of Ref. [24] . 

The 31 SDs of the a-cluster wave functions are used 
in the GCM calculation in Ref. [24]. We place the a- 
particle wave functions at the same positions as those of 
Ref. [24]. In Ref. [24], the single-particle wave functions 
of the SDs are the Gaussian wave packets. Instead of the 
Gaussian wave packet, we employ the HF orbitals of the 
a particles. In Fig. 11, we show density distributions of 
selected SDs among those 31 SDs. 

In Fig. 12, we show excitation spectra from configura- 
tion mixing calculations using the 31 SDs. The left panel 
shows our calculation using SLy4 interaction. The mid- 
dle panel shows the GCM calculation using Volkov No. 
I force [24]. The results for the ground rotational band 
are similar to each other. In fact, in both calculations, 
the moment of inertia is too large. The O2 state appears 
at around 7 MeV in our calculation, slightly lower than 
that of Ref. [24]. 

In the parentheses in Fig. 12, we show the calculated 
binding energies in the ground state. The absolute val- 
ues of the binding energies is very different between our 
calculation and that of Ref. [24] . In our calculation using 
SLy4 interaction, the binding energy is -75.1 MeV and is 
much smaller than the value shown in Fig. 2. A major 
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FIG. 11. Density distributions of some SDs out of 31 SDs which are used in Ref. [24]. Unit of vertical and horizontal axes are 
fm. 



part of the difference comes from the spin-orbit interac- 
tion which contributes little in the calculation using the 
alpha-cluster wave functions only. 

We next perform a configuration mixing calculation 
employing both the 50 SDs prepared by the imaginary- 
time method and the 31 SDs of three-a configuration. 
In Fig. 13, we compare the three calculations: the con- 
figuration mixing calculation using 50 SDs prepared by 
the imaginary time method (IT), the configuration mix- 
ing calculation using 31 SDs of three-a configuration 
(3a), and the configuration mixing calculation using both 
(IT+3a). The calculation labeled by three-a is the same 
as that shown in the left part of Fig. 13, except that the 
total energies are plotted here. 

After mixing both configurations of the imaginary-time 
and the three-a, we find the results are very close to the 
calculation using the imaginary-time configurations only. 
Namely, 31 SDs of the three-a wave functions do not mix 
with those prepared by the imaginary-time method. This 
is due to the large energy difference between those two 
sets of configurations. 

In the calculation using configurations generated by 
the imaginary-time method, the contribution of the spin- 
orbit interaction to the binding energy is as large as 17 
MeV with SLy4. This large energy gain is missing in the 
pure three-a configurations. 

In Table VII, we show the calculated radii and the 
M{EO) transition strength. Using the 31 SDs of three- 
a wave functions, our calculation gives large values for 
both the 0^ and Oj states. The radius of the Oj state 
is 3.31 fm, close to the value by the GCM calculation, 
3.4 fm. The M{EO) value is also large, 8.72efm , even 
larger than the three-a GCM calculation [24]. However, 
in the configuration mixing calculation using both con- 
figurations, our calculated values are very close to the 
calculation using the 50 SDs prepared by the imaginary- 
time method. This result is consistent with the fact that 
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FIG. 12. Energy levels of ^^C employing 3a SDs with Skyrme 
SLy4 interaction (left), the results of GCM calculation of 
Ref. [24] (center) and experiments (right). 



the energy spectra is very little affected by adding the 
three-a wave functions. 



V. SUMMARY 

We have investigated structure of the ^^C nucleus em- 
ploying a configuration-mixing approach with Skyrme in- 
teraction. In this approach, we first generate a number 
of Slater determinants using the imaginary-time method 
[19] starting from initial Slater determinants prepared 
in a stochastic way. These Slater determinants show 
various shapes and clustering. They are projected on 
parity and angular momentum, then, are superposed 
by the configuration-mixing calculation. We have gen- 
erated several sets of Slater determinants and compare 
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EXP 


IT 


IT + 3a 3a 


3a(Uegaki) 


radius(O^) 

radius(Oj) 

M{EQ;0+ ^0+) 


2.31±0.02 
5.4± 0.2 


2.53 
2.76 

4.57 


2.54 2.80 
2.73 3.31 
4.13 8.72 


2.40 
3.40 
6.6 



TABLE VII. Radii (fm) and M{EO) (efm^) calculated in various model spaces. Results of GCM calculation [24] is also shown. 
See text for details. 
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FIG. 13. Energy levels of ^^C in the configuration mixing cal- 
culation with the SDs constructed by the imaginary-time evo- 
lution (IT), and 3a, and all of these configurations (IT-|-3a). 
See text for details. 



tures of the structure of ^^C. The energies and the B{E2) 
transition strength of the ground state rotational band 
are well described. The lowest excited states of negative 
parity, 1~, 2~, and 3~, are also reasonably described, 
although the excitation energies are slightly too high. 
The Slater determinants which dominate in these states 
show three-a structure. Our calculation also reproduces 
the excitation energy of the Hoyle (Oj) state reasonably. 
This state is found to be described by superposition of 
many Slater determinants, consistent with the former 
cluster-model calculations. However, the radius of the 
0;^ state is calculated to be significantly smaller than 
those. The energy gain associated with the spin-orbit in- 
teraction in the present method seems to be responsible 
for this difference. The three-a liner-chain structure ap- 
pears at around 15 MeV excitation, forming a rotational 
band. 

The success for the description of ^^C nucleus reported 
in this paper clearly shows that the present approach 
is promising for systematic description of various many- 
body correlations including clustering in light nuclei. 



the results with the different sets, to quantify the re- 
liability of the calculation. A few low-lying states for 
each parity and angular momentum are well converged 
with small variance among the different sets of the Slater 
determinants. This fact indicates that the present cal- 
culation provides unique and convergent results for the 
ground and a few low-lying excited states, once the effec- 
tive nucleon-nucleon force, the Skyrme interaction in the 
present calculation, is given. 

Our calculation reasonably reproduces the overall fea- 
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